OG 2.1.03 

Hour-Scale Multi wavelength Variability in Markarian 421 

M.Catanese 1 , LH.Bond 2 , S.M.Bradbury 2 , A.C.Breslin 3 , J.H.Buckley 4 , A.M. Burdett 2 5 , 
D.A. Carter-Lewis 1 , M.F. Cawley 6 , S. Dunlea 3 , M. D'Vali 2 , D.J. Fegan 3 , S.J. Fegan 5 , J.P. Finley 7 , 
J.A. Gaidos 7 , T.A. Hall 7 , A.M. Hillas 2 , D. Horan 3 , J. Knapp 2 , F. Krennrich 1 , S. Le Bohec 1 , 
R.W. Lessard 7 , C. Masterson 3 , P. Moriarty 8 , J. Quinn 3 , H.J. Rose 2 , F.W. Samuelson 1 , G.H. Sembroski 7 , 
V.V. Vassiliev 5 , T.C. Weekes 5 , L. Maraschi 9 , G. Fossati 10 , F. Tavecchio 9 , L. Chiappetti 11 , A. Celotti 10 , 
G. Ghisellini 9 , P. Grandi 12 , E. Pian 13 , G. Tagliaferri 9 , A. Treves 14 
1 Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, USA 
^Department of Physics, Leeds University, Leeds, LS2 9JT, UK 
3 Experimental Physics Department, University College, Belfield, Dublin 4, Ireland 
^Department of Physics, Washington University, St. Louis, MO 63130, USA 
5 Fred Lawrence Whipple Observatory, Harvard-Smithsonian CfA, Amado, AZ 85645, USA 
qv | 6 Physics Department, St. Patrick's College, Maynooth, County Kildare, Ireland 

t— l ' 7 Department of Physics, Purdue University, West Lafayette, IN 47907, USA 

£^ ■ ^Department of Physical Sciences, Galway-Mayo IT, Ireland 

9 Osservatorio Astronomico di Brera, Milano, Italy 

W SISSA/ISAS, Trieste, Italy 

ll IFC/CNR, Milano, Italy 

12 IAS/CNA, Roma, Italy 

13 TESRE/CNR, Bologna, Italy 

^^ ■ 14 Universita dell 'Insubria, Como, Italy 

o: 

£J ' Abstract 

o: 

Q\ , Markarian 421 was observed for about four days with BeppoSAX and the Whipple Observatory 7-ray tele- 

scope in April 1998. A pronounced, well-defined, flare with hour-scale variability was observed simultane- 

S^. ■ ously in X-rays and very high energy 7-rays. These data provide the first evidence that the X-ray and TeV 

intensities are well correlated on time-scales of hours. While the rise of the flare occurred on a similar time- 
scale in the two wavebands, the decay of the flare was much more rapid in 7-rays, providing the first clear 
indication that the X-ray and 7-ray emission may not be completely correlated in Markarian 421. 

1 Introduction: 

Blazars are a class of active galactic nuclei whose emission is believed to arise predominantly from a 
relativistic jet whose axis makes a small angle with our line of sight. More than 60 blazars have been detected 
with the Energetic Gamma-Ray Experiment (EGRET) (Hartman et al. 1999). A few BL Lacertae objects (BL 
Lacs), a sub-class of blazars, have also been detected as TeV 7-ray emitters (Ong 1998). No model for the 
origin of the 7-ray emission is generally accepted at this time. Two popular classes are those in which high 
energy electrons produce the 7-rays by inverse Compton scattering of low energy photons (e.g., Maraschi, 
Ghisellini & Celotti 1992; Dermer, Schlickeiser & Mastichiadis 1992; Sikora, Begelman & Rees 1994) and 
those in which high energy protons produce 7-rays by initiating cascades in the jets (e.g., Mannheim 1993). 
Contemporaneous observations at several wavelengths can be used to derive physical conditions in and around 
the blazar jet and may resolve which emission mechanism operates in the objects. 

Markarian 421 (Mrk421) is the closest known BL Lac (z = 0.031) and is an established very high energy 
(VHE, E>250GeV) 7-ray source (e.g., Punch et al. 1992; Petry et al. 1996). Mrk421 is also an EGRET 
source (Hartman et al. 1999). The VHE emission from Mrk421 is extremely variable, showing flaring activity 
on time scales as short as 15 minutes (Gaidos et al. 1996) with little or no baseline level emission (Buckley 
et al. 1996). The spectrum of Mrk421 is consistent with a power law that extends to at least 10 TeV with no 
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evidence of a sharp cut-off, and no evidence of variability (e.g., Krennrich et al. 1999a,b). Multiwavelength 
campaigns on Mrk421 (e.g., Buckley et al. 1996) have revealed correlations between X-rays and VHE 7-rays 
and evidence for correlated optical/UV variability. The flux amplitude of the X-ray and VHE 7-ray variations 
was similar and the variability time profiles were the same, on day-scales. These rapid, correlated variations 
have permitted stringent limits to be placed on the Doppler factor and magnetic fields of the Mrk421 jet (e.g., 
Buckley et al. 1996) and these data have become important tests for emission models (e.g., Mannheim 1998; 
Buckley 1998; Tavecchio, Maraschi & Ghisellini 1998). 

Despite the successes of these campaigns, the light curves were not densely sampled, so the multiwave- 
length variability could not be measured on time-scales less than one day. In order to better measure the flaring 
behavior of Mrk421, several more intense multiwavelength campaigns were conducted in 1998 using longer 
exposures in X-rays and VHE 7-rays, and combining the data from several VHE 7-ray telescopes. Here, we 
present the results of a campaign in 1998 April with BeppoSAX and the Whipple 7-ray telescope. 



2 Observations: 

2.1 BeppoSAX: The scientific payload carried by 
BeppoSAX is fully described in Boella et al. (1997a). 
The data of interest here derive from three coaligned in- 
struments, the Low Energy Concentrator Spectrometer 
(LECS, 0.1-10 keV, Parmar et al. 1997), the Medium En- 
ergy Concentrator Spectrometer (MECS, 2-10 keV, Boella 
et al. 1997b) and the Phoswich Detector System (PDS, 12- 
300 keV, Frontera et al. 1997). 

The observations with BeppoSAX reported here con- 
sist of two exposures lasting approximately 100 kilosec- 
onds each. The data reduction for the PDS was done using 
the XAS software (Chiappetti & Dal Fiume 1997), while 
for the LECS and MECS linearized cleaned event files 
generated at the BeppoSAX Science Data Centre (SDC) 
were used. No appreciable difference was found extract- 
ing the MECS data with the XAS software. Light curves 
were accumulated from each instrument with the usual 
choices for extraction radius and background subtraction 
as described in Chiappetti et al. (1999). 

2.2 Whipple: The VHE 7-ray observations were 
made with the Whipple Observatory 10 m telescope (Caw- 
ley et al. 1990). At the time of these observations, the 
telescope camera consisted of 33 1 photomultiplier tubes 
with a combined field of view of 4.8°. Also, light-cones 
were not in place and this, as well as reduced reflectiv- 
ity of the mirrors, resulted in a somewhat higher energy 
threshold than usual for the telescope, ~500 GeV Events 
were parameterized with a standard moment analysis and 
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Figure 1 : The 7-ray and X-ray light curves of Mrk 
421 from 1998 April 21 to 24 (UT). The top panel 
shows data for E>2TeV The middle and bottom 
panels show BeppoSAX data for the three instru- 
ments with energy ranges specified in the plots. 
All of the count rates are normalized to their re- 
spective averages (listed in the panels) during this 
observation period. 



candidate 7-rays were selected using a variation of the Su 

percuts analysis (Reynolds et al. 1993) appropriate for the large camera field of view and for maintaining a 

constant energy threshold as a function of observation elevation (see Tableland discussion below). 

Observations were taken on the nights of 1998 April 21, 22, 23 and 24. To permit longer observations 
within each night, data were taken over a large range of zenith angles (~7° to 60°). The collection area and 



Table 1 : Elevation dependent parameter cuts for Mrk 42 1 analysis 



Parameter 


<35° 


Observation zenith angle 
40° 45° 


55° 


a 

length 

width 

distance 

asymmetry 

length/size 

max2 


<7.5° 

> 0.10°, < 0.45° 

> 0.04°, < 0.18° 

> 0.30°, < 1.60° 

> 0.0 

>80d.c. 


<7.5° 

> 0.10°, < 0.45° 

> 0.04°, < 0.15° 

> 0.30°, < 1.60° 

>0.0 

> 80 d.c. 


<7.5° 

> 0.10°, < 0.38° 

> 0.04°, < 0.13° 

> 0.30°, < 1.60° 

> 0.0 

> 80 d.c. 


< 10.0° 

> 0.10°, < 0.30° 

> 0.04°, < 0.12° 

> 0.30°, < 1.60° 

>0.0 

< 0.00078°/d.c. 6 

> 60 d.c. 


size 
Effective area ratio d 


> 1300 d.c. 
1.0 


> 900 d.c. 
1.5 


> 600 d.c. 
2.6 


> 200 d.c. 
4.7 



a \ d.c. = 1 digital count. 

6 Size cut required for energy threshold of 2 TeV. 

c Ratio is relative to the effective area at a zenith angle of 20° 



energy threshold increase with zenith angle and the 7-ray selection is a function of zenith angle, so the ob- 
served 7-ray rates can change as a function of elevation even from a source of constant 7-ray emission. To 
obtain a light-curve which shows only the intrinsic source variations, it was necessary to determine software 
cuts which result in a constant energy threshold as a function of elevation and to calculate collection areas as 
a function of elevation for those energy thresholds in order to normalize the 7-ray rates. Because contempo- 
raneous observations of the Crab Nebula were not available with sufficient statistics over all the zenith angles, 
the analysis presented here relies entirely on Monte Carlo shower simulations. 

Simulated 7-ray induced showers at zenith angles of 20°, 40°, 45°, and 55° were generated with ISUSIM 
(Mohanty et al. 1998) to determine the size cut required at each zenith angle to obtain a common energy 
threshold of 2 TeV and to estimate the collection areas at these zenith angles for the 2 TeV energy threshold. 
To normalize the >2TeV rate measurements at the different elevations we multiply them by the ratio of 
effective collection area at a given zenith angle to the effective collection area at a zenith angle of 20°. The 
software trigger threshold applied at each zenith angle to set the energy threshold at 2 TeV and the ratio of 
effective areas for the four zenith angle ranges are shown in Table Q The results reported here are based on 
limited statistics and are therefore preliminary. The aim of this analysis is to derive normalized fluxes as a 
function of time rather than absolute fluxes and energy spectra. 

3 Results: 

The light curves for the 7-ray and X-ray observations are shown in Figure^ Three X-ray energy bands are 
shown and the 7-ray light curve shows the normalized E>2 TeV data for the measurements. Each 7-ray point 
represents a 28 minute observation. The count rates for the measurements are normalized to their respective 
averages for these observations. The rise and fall of a large amplitude flare is clearly evident in all data sets 
on the first day of observations. Observations after the first day did have detectable fluxes, but showed no 
significant variability on day or shorter time-scales. For the observations on April 21, the amplitude of the 
X-ray flaring increases with increasing energy, but it is close to a factor of 2 in all three bands. The >2TeV 
light curve shows a 4-fold variation in flux. The flux in the 0.1-2keV, 4-10keV, and >2TeV energy bands 
peaks at approximately the same time (within one-half of one hour), but the decay time for the TeV light curve 
is significantly shorter than that of the LECS and MECS light curves. The 12-26 keV light curve measured 
by the PDS instrument appears to peak slightly later than the others, but the statistical uncertainty in the data 



precludes a definitive measurement. A detailed investigation of possible leads or lags in the data is underway. 
The TeV spectrum does not change significantly during the rapid flare, nor is it significantly different than 
previous measurements (Krennrich et al. 1999b). 

4 Discussion and Conclusions: 

These observations show, for the first time, that the TeV and keV fluxes from Mrk421 are correlated on 
hour time scales while at the same time indicate that the 7-rays and X-rays are not completely correlated. 
Neither the larger variability amplitude at TeV energies than at X-rays nor the difference in the variability time 
scales have been seen previously. The reason for the difference in the decay time-scale of the flare at TeV 
and keV energies is not clear. The differences could reflect the nature of the flaring mechanism. For example, 
a variation in the electron spectrum and the energy density of the low energy photons up-scattered to TeV 
energies (Maraschi et al. 1999) might produce such a flare. The differences may also indicate that the particles 
which produce the X-rays are not the same as the particles which produce the TeV 7-rays. This is possible 
in models where the progenitor particles are electrons or protons. In addition, the region of the broadband 
spectrum of Mrk421 observed by BeppoSAX spans the end of the synchrotron emission and the onset of the 
high energy emission (c.f., Buckley et al. 1996). As such, the X-ray emission may reflect contributions from 
more than one population of source particles, regardless of the emission mechanism. Detailed model fitting of 
this data, which is beyond the scope of this paper, is necessary to investigate these possibilities. 
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